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Abstract 

Zee model leads naturally to two heavy, strongly degenerate and almost 
maximally mixed neutrinos and to one light neutrino with small mixing. This 
pattern coincides with the one needed for a solution of the atmospheric neu- 
trino problem by i/„ — * u T oscillations and for existence of the two component 
hot dark matter in the Universe. Furthermore, the oscillations Va — > V e can 
be in the range of sensitivity of KARMEN, LSND experiments. Phenomenol- 
ogy of this scenario is considered and possibility to check it in the forthcoming 
experiments is discussed. Scenario implies large values and inverse flavour hier- 
archy of the couplings of the Zee boson with fermions: f er <C f^ T < f e ^ ~ 0.1. 
Main signatures of scenario are: strongly suppressed signal of f« — > v T oscilla- 
tion in CHORUS and NOMAD experiments, so that positive result from these 
experiments will rule out the scenario; possibility of observation of v e — > v T os- 
cillations by CHORUS and NOMAD; corrections to the muon decay, neutrino- 
electron scattering at the level of the experimental errors; branching ratio 
B(u — > e~f) bigger than 10~ 13 . The solar neutrino problem can be solved 
by introduction of additional very light singlet fermion without appreciable 
changes of the active neutrino pattern. 
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1 Introduction 



Zee model [0 is the simplest model which explains a smallness of neutrino masses by 
physics at the electroweak scale. It can be considered as an alternative of the see-saw 
mechanism 0. 

Zee model suggests the existence of charged scalar field h, singlet of the SU(2), 
and two doublets of the Higgs bosons $1, $2- The singlet h couples to lepton doublets 
^il — (vi, l~), (I = e, //, r) as well as to Higgs doublets; leptons are assumed to couple 
to doublet $i only. The appropriate terms in the Lagrangian are 

Czee = fw^L^t'Lh + C 12 <$>1 iT 2 <$> 2 h ] + -^r%^ll R + h.C, (1) 

where C\ 2 = —c 2 \ are real mass parameters, the couplings fu> are antisymmetric in I 
and £'. The interactions ([!]) generate neutrino mass terms in one loop. 

Zee model gives very distinctive pattern of neutrino masses and mixing. For not 
too strong hierarchy of the couplings fty the two heavy neutrinos, Vi, ^3, are strongly 
degenerate and mix almost maximally in u^, u T . The first neutrino v\ practically 
coincides with v e and has much smaller mass: 

mi m 2 ~ 777.3. (2) 

It was marked |§ that this pattern coincides with the one needed to solve simultane- 
ously the atmospheric neutrino problem || and the problem of the hot dark matter 
in the Universe ||. Indeed, the deficit of the atmospheric muon neutrinos can be 
explained by the oscillations — u T with practically maximal mixing. Two heavy 
neutrinos with masses vn 2 ~ 777,3 ~ (1 — 5) eV compose two component hot dark mat- 
ter (which may give even better fit of the cosmological data than one component) ||. 
Furthermore, the oscillations — > v e and z7 M — > V e can be at the level of sensitivity 
of existing experiments: BNL 0, KARMEN (see 0). Later it was marked g [|] 
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that the model can immediately accommodate positive LSND result [|K|. 

In this paper we will consider phenomenology of the outlined scenario, and in 
particular, the possibility to check it by forthcoming experiments. In sect. 2 we 
describe the scenario in details. Sect. 3 is devoted to oscillations. In sect. 4 we find 
the bounds on the Zee coupling constants. In sect. 5 implications of data on the 
muon decay, neutrino electron scattering, e — /i — r universality to the scenario are 
considered. Predictions for \i — > e^y and z/ 3 ( 2 ) — > v^y are given. In sect. 6, we describe 
a modification of the Zee model which is able to solve the solar neutrino problem. 
Sect. 7 contains our conclusions. 

2 Scenario 

The neutrino mass matrix of the Zee model in flavor basis, v = (z/ e , z/ r ), can be 
written as 

/ e sin#\ 
M = m \ e cosfl , (3) 
ysin^ cos 6* J 

where mo is the basic mass scale. Mixing angle 9 and parameter e can be naturally 
much smaller than 1. (We will discuss the relation of these parameters with the 
parameters of the Lagrangian (p]) in sect. 4.) 

In the case cos 9 ^> sin 9, e the eigenvalues of matrix ([|) are 

mi = —m Q esm29 , m 2 ,3 = m (±l — -esin2#) , (4) 

and the mixing matrix S which diagonahzes (|3|) is 

j / \pl cos 9 sin 9 + e cos 9 sin 9 — e cos 9 \ 
S~—\-^2sm9 cos9 cos9 , (5) 

^{ -V2e 1 -1 ) 
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(uf = Si/, where v = {v\, z/ 2 , z/ 3 ) are the mass eigenstates). According to ([4[) the 
states V2 and z/3 are approximately degenerate, and their masses (~ mo) are much 
larger than the mass of v\. The mass squared difference is 

Am^ = 2e sin 2Qm\ <C m\ , (6) 

where Amfj = m 2 — m|. For z/i component model gives Am^ ~ Am^ ~ Wq and the 
ratio of mass differences equals 

^k = 2esin29. (7) 
Am 2 21 

Thus, the Zee mass matrix gives two different scales for the mass squared differences 
and the maximal mixing between two heaviest neutrinos. 

As it was outlined in the introduction we will consider the following scenario: 

• V2 and z/3 form the two component hot dark matter so that 

m = m HDM = (f - 5) eV. (8) 

• The z/ M — z/ T oscillations with practically maximal depth explain the atmospheric 
neutrino deficit and therefore 

Am 2 2 = Am 2 atm ~ (0.3 - 3) • fCT 2 eV 2 . (9) 

• The oscillations — v e and — V e with Am\ x « m< HDM can be in the re- 
gion of sensitivity of the KARMEN and LSND experiments. For ttlq in the 
cosmologically interesting domain (|8]) this means that 

sin 2 26 < sin 2 26 eiM ~ (1 - 3) • f (T 3 , (f 0) 

where 9 e ^ is the experimental bound (or preferable value in the case of positive 
result) for the u e — z/ M mixing angle. 



Substituting m and sin 2 29 from (||) and ( [10] ) in (|6]) we get 

According to (|) and QTT|) e = 1(T 3 - 0.5, with typical value 3 • 10~ 2 . Thus all 
oscillation parameters of the model (mass squared differences and mixing angles) can 
be fixed by the experimental data immediately. 



3 Neutrino Oscillations 

In terms of the elements of the mixing matrix (|5]), S a i, the oscillation probability can 
be written as 

P[v a -> vp) = 5 a p - 4 SaiSpiSajSpj sin 2 ^ ^ , (12) 

where E v is the neutrino energy and L is the distance. We neglect the CP violation, 
so that the elements S ai are real. Let us consider the probabilities (0) for short and 
long distances separately. 

1). In the short distance limit the phase difference due to Am 2 2 is small: A ^ 2jL <C 
1. Taking into account that Am| x ~ Am^ and using matrix (|5|) we find 

-> v e ) w P(V^ -> 77 e ) = sin 2 29 sin 2 (^|^J • (13) 

This result is applied to E776 [§, KARMEN and LSND [|Kj. 

For Ufj, — v T we get 

P(„ _> „ T ) = 4 ( e sin fl) 2 sin 2 f ^f^l + cos 2 sin 2 f . (14) 

The ^ — z/ r oscillations with large mass splitting are doubly suppressed because of 
sin 9 <C 1 and e <C 1. This smallness is related to the fact that according to @ the 
admixture of v\ in v T is suppressed by e. The mode of oscillations with the smallest 



mass splitting (second term in (fH|)) may give a comparable contribution. For values 
of sin 2 29 and e from (|TD|) and ( |TTD correspondingly we obtain 

P(v„ -> i/ T ) ~ 10- 7 - 10- 5 . (15) 

If both e and Am 2 are near the upper bounds the probability can be as big as 10~ 4 . 



These values are still below the sensitivity of CHORUS and NOMAD fill , but they 



may be in the regions of sensitivity of planning experiments E803 at Fermilab and 



E889 at BNL 12 



For v e — v T channel we get 



P(z, e ^z, T )=4( e cos^) 2 sin 2 [^^') . (16) 



If e > 10 1 and m > 4 eV, then P(u e — > v T ) ~ 10 2 - 10 1 and the v e — > i/ T 
oscillations can be detected by CHORUS and NOMAD. 

Thus the observation of signals of the v e — > v T oscillation and absence of signal 
from — > v T mode in CHORUS and NOMAD are the signatures of the Zee model. 
The scenario under consideration will be ruled out if CHORUS and NOMAD find 
signals of the —>■ v T oscillations. 

According to @ for fixed Am 2 2 and mo, the parameter e is inversely proportional 
to sin#. Therefore P[y e — > u T ) increases when Piv^ — > v e ) decreases, as is shown in 
Fig.l. 



Fig. 1 



In particular, if P(v^ — > V e ) ~ 3 • 10 -3 (the level of the LSND result) and m 2 , > 6 
eV 2 , then P{y e — > v T ) < 3 • 10 -5 which is beyond the sensitivity of the CHORUS and 
NOMAD. On the contrary, for P{y e — ► z^ T ) > 2 • 10~ 2 which can be observed by these 
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experiments one has P[y^ — > v e ) < 10 5 . Thus a comparison of results from searches 
for P(v e — > z/ T ) and P[y e — > z/ M ) oscillations can give crucial check of the model. 

For the parameters under consideration there are strong resonance transitions 
v e — > Vp, v T and — > u e in the inner parts of the collapsing stars. As the consequences 
one predicts: (i) disappearance of the neutronization peak, (ii) hard v e spectrum at 
the cooling stage, (hi) additional energy release in the inner parts of star which 
will stimulate shock wave revival desired for the star explosion, (iv) In the same 
time the — > u e conversion leads to suppression of the r-processes responsible for 
nucleosynthesis of heavy elements unless itlq < 2eV ||13| . 

2). In the long distance limit experiments are sensitive to oscillations stipulated 
by small mass difference Am^ and the oscillations due to large mass difference are 
averaged out. We get the results 

P{v.^p e ) = i(sin2e) 2 -cos 2 ^(sin 2 ^-e 2 cos 2 ^)sin 2 (^i^) , (17) 

2 4:E U 

P{y^v T ) = 2( e sin^) 2 + cos 2 ^sin 2 ^^ 2 ^^ , (18) 

which are applied to the atmospheric neutrinos. Notice that P{vn — > v e ) is suppressed 
due to both sin# <C 1 and e C 1, and the dominant effect comes from — > v T 
oscillations as we suggested in the introduction. 



4 Parameters of Zee Model 

In terms of parameters of the Lagrangian (P the elements of the mass matrix (H) 
equal 

tanfl = ^, 6= . /eM (Hhi)\ (19) 
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»o^;V^^=^S- (20) 



and |3 

-g sin 20 cot /3 Mf 
64 V / 2M W 7T 2 n Mf 

Here m r is the tau lepton mass, g is the weak coupling, my/ is the H^-boson mass, 
tan j3 =< $i >o / < $2 >o is the ratio of the VEV of two Higgs doublets. The angle 
j3 determines physical charged Higgs boson: $ + = $^ cos (3 — sin /3, where Qf, $2" 
are two charged Higgs fields from the doublets. The angle <ft is the mixing angle of 
the Zee singlet and the physical charged component of the Higgs doublet $ + : 

h = cos <pHi + sin <pH 2 , 
$+ = - sin (j)H x + cos (j)H 2 , (21) 

where H\ and i?2 are the eigenstates of the mass matrix with masses M\ and M 2 , 
and the mixing angle is determined by 

tan 24, = Wfr-'<*.Jfr (22) 

y'fM? - M'i f - (4 v / 2 9 - 1 ci2M„,) 2 

As we have seen in sect. 2, the parameters of the mass matrix (|3]) m , e, 9 can 
be fixed by the data. This in turn allows one to determine the ratios of the constants 



fij using (ID 



^ = tan^<l, (23) 



and 



fe/i _ Am ltm ( m r 



. : . (24) 

For sin 2 29^ = 2 x 1(T 3 , Am 2 HDM = 6 eV 2 and Am 2 tm = 10" 2 eV 2 Eq. @ gives 
fe/i/f/iT = 5.3 which means an inverse hierarchy of the Zee boson couplings with 
f eil being the largest one ||. For fixed value Piy^ — > v e ) the mixing angle 9 is the 
function of Am 2 x = m\. Using this dependence we get from fl2"3|) and (E3j) the ratios 
fefj,/ fer and f^r/ fer as the functions of mo for fixed value of -P(^ M — > f e ) (see Fig. 2). 
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For Piv^ — > z/ e ) = 1.5 x 10~ 3 , (which is in the range of sensitivity of KARMEN and 
LSND) we find f efl ~ / MT ^> / er at mo = 5 eV. This relation may testify for certain 
horizontal symmetry. Below m = 5 eV, there is an inverse flavour hierarchy of the 
couplings, f etl > f^ T ^> f er . For Piy^ — > u e ) < 10~ 4 one gets the inverse flavor 
hierarchy already below m = 10 eV. 



Fig. 2 



The absolute value of the coupling constants can be fixed by (pOP- For values of 
parameters: sin0 ~ O(10 _1 ), tan/3 ~ O(10), Mi ~ M 2 ~ O(500GeV) we get: 
/ eM = 10 -2 — 1. That is the scenario implies quite big couplings constants of the Zee 
boson. 

5 Constraints from the Electroweak Processes 

Since the constants f efl , / MT are rather big the Zee singlet can give observable contri- 
butions to different weak processes. The effective four-fermion Lagrangian induced 
by the Zee boson exchange can be written (after appropriate Fiertz transformation) 

as 

% [*V/*(1 - j 5 )eV^(l - 75 )e - u^(l - 75 )/^e7 M (l - 7s)e + ...] , (25) 



where 

=2 | , (26) 



and 

1 cos 2 6 sin 2 



1 f 2 
V2G F W 
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Notice that only usual left handed components of leptons participate in the interac- 
tions with Zee boson, and therefore the Lagrangian (|2"5| ) has usual V — A form. 

In the case neutrino electron scatterings, e~ — > V ^ e~ and e~ — > e~ , the 
contribution from ( pop leads to a change of the g e L coupling: g e L — > g e L + £ . CHARM 



II experimental data on gi and gn |T5[ agree well with predictions of the Standard 
Model. Therefore £ should be smaller than the experimental error Ag£: £ < Ag£. 
Using (|26|) we have explicitly 

^ < 0.036GV . (28) 
M 

The Zee singlet exchange leads also to the lepton number violating process e~ — > 
v T e~ which contributes to v^e scattering incoherently. Its amplitude is proportional 

t° Jenifer- 

The second term in the Lagrangian ( |25| ) gives the renormalization of the four 
fermion coupling Gp of the muon decay. Assuming that the effect of the Zee boson 
on the decay rate is smaller than 0.1%, (so that it does not destroy the agreement in 
the electroweak precision tests) we find 

^<7-l(T 4 G F . (29) 
M 

Also the modes of the muon decay with lepton number violation appear: fx — > v T e V ei 
fi — > e I7 T , fi — > v T e I7 M which contribute to the total decay rate incoherently. 

The result (^) allows one to get the bounds on masses and mixing of scalar 
bosons. Indeed, using expression for the mass (p0|) we can find f efl as the function 
of (f), P and Mi. Substituting f efl = f efJi ((fi, j3, Mi) into (p9j) we find the lower bound 
on sin0 as the function of Mi for different values of Piv^ —>■ v e )i tan/5 and 
fixed M2 = 300 GeV (see Fig. 3 (a) — (c)). Notice that the most strong bound is for 
Mi = Mi- Forbidden region becomes larger with increase of tan/3 as well as with 
decrease of itlq and P{p i1 —>■ v e )- Big region of parameters exists in which all the 
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restrictions are satisfied. Furthermore, one of the charged Higgses can be at the level 
of lower kinematical bound. 

~~Fig. 3 (a), 3 (b), 3 (c) 

The bound on the model follows also from e — \x — r universality. The ex- 
pected deviation from universality due to the Zee boson contribution is |1 — g T j 1 g^\ ~ 
(fe )/(GpM ), where and g T are the weak coupling constants of the charged 
currents with [i and r. Recent measurement of the branching ratio of the decay 
r — > ev e v T at OPAL [[U| gives the ratio of couplings g T / g^ = 1.0025 ± 0.0060, and 
the corresponding bound on the parameters of model is weaker than (f29|). 



The model leads to the radiative decays of the muon /i — > and neutrino 
^3(2) - * through the one-loop diagram with Zee singlet. 

The branching ratio of the /i — > is 

B( ^^<^)(wtJ ■ (30) 

Using (p!9|), p0| ) and (|27D we can express it as B(fi — > ej) = A(sm4>, Mi, tan/3)r?7.Q. 
The branching ratio becomes smaller with increase of sin and decrease of tan j3 (see 
fig. 4). The present experimental upper bound B < 4.9 x 10 -11 |l7j (shown by the 
horizontal dashed line) will be strengthen soon up to 5 x 10~ 13 by the experiment 
at MEGA in LAMPF(Los Alamos). Future experiment \T$\ will push the limit to 
3 x 10~ 14 . The results from these experiments combined with bounds from precision 
tests (fig. 3) will cover essential part of the parameter space of the model. 



Fig. 4 
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The life time of Vi — > 1/17 (z = 2, 3) equals U [|14 

^ 2 ' a 



TV; 



am,- 



m 



2— £C, ; 1 - cos 20 



77? 



2 w ^ 



a, 



1 rai 

777-i 



(31) 



where 



In 



M2 H2 



In 



777? 



(2 



/i,r 



(32) 



The life time r(ui —>■ v{) depends mainly on the charged Higgs scalar masses Mi and 
M2; fit' and sin enter only via the mass of neutrino. For mo = 1 — 10 eV the life time 
is in the interval 10 22 — 10 29 years. This may have some cosmological implications. 

In the limit of f eT = the anomalous magnetic moment of neutrino which corre- 
sponds to fl3lD equals Jl4 

fMu ~ — AemoCr , (33) 

where C T is defined in fl52"p. For Mi ~ M 2 ~ 300 GeV and m = 2.65 eV, we get 
fi u ~G x 10~ 16 e/2m e . 



6 Solar neutrinos 



For solar neutrinos all oscillations are averaged and from ([□) one gets survival prob- 
ability 

P(y e _> !/ e ) = cos 4 9 + ~ sin 4 + 0(e 2 ) . (34) 

There is no dependence of suppression of the v e - flux on energy and for e, sin 2 9 <C 1 
the effect is small. Thus in the considered scenario there is no solution of the solar 
neutrino problem. 

Let us suggest that apart from three known neutrinos also singlet (because of 
the LEP bound) neutrino v s exists. This neutrino mixes with electron neutrino so 
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that the resonance conversion u e — > v s explains the deficit of the solar z/ e -flux. The 
explanation requires the mass squared difference and the mixing angle in the intervals 

Am 2 = (4 - 10) • l(T 6 eV 2 , sin 2 29 es = 1(T 3 - 10~ 2 . (35) 

The singlet neutrino could be the right handed counterpart of the known neutrino 
components or new very light fermion which comes from some other sector of theory. 

The mass of the lightest neutrino in the Zee model (which is essentially the v e ) is 

A 777^ 

mi = m e sin 29 « ^ ~ (1 - 5) ■ 10~ 3 eV . (36) 

2rriHDM 

Squared mass m\ is close to Am 2 desired for solar neutrinos (|35|) . This means that 

the mass of singlet neutrino, m s , should be rather close to m% (recall that for the 

resonance conversion one needs m s > mi): 

m a - mi _ Am 2 
mi 2m 2 

For mi > 4 • 10~ 3 eV one gets from this equation Am/ mi < 0.2 . 

Let us consider the simplest scheme with only one singlet neutrino. We extend 
the Lagrangian of the Zee model by adding the terms: 

0i®v s + m ss v T s v s . (38) 

All couplings fi can be of the same order. The first term leads to mixing of the v s 
with the active neutrinos: m; s = fi(&). Performing block diagonalization of 4 x 4 
mass matrix we get the mass matrix for the {y s — v e ) system: 

M _ / m ss m es \ 
\ m es mi J 

where m es = f e (&) arid mi is fixed in fl55|) . The mixing angle is then 

sin2g es « 2mes . (40) 
m ss - mi 
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If the mass m ss is not too close to mi, we get 

m s ~ (2 - 4) • l(T 3 eV , m es ~ l(T 4 eV. (41) 
With increase of m s (and consequently the degeneracy) m es can be further diminished. 

According to flilD a solution of the solar neutrino problem implies very small 
Yukawa coupling f e < 1CT 15 which is of the order m EW /m string , where m string ~ 10 18 
GeV is the superstring scale. The mass of the singlet neutrino m s is of the order 
m 3/2/ ' m string which also may indicate the SUSY origin of the singlet. 

Mixing of the singlet neutrino with high mass states z/ M , v T is of the order 

sin 2 26 ns ~ sin 2 26 es ■ sin 2 9 ~ 1(T 7 , (42) 
so that the bound from the primordial nucleosynthesis [^TJ can be satisfied. 

The influence of the singlet fermion on "standard" structure of the Zee model is 
negligibly small and the results of the previous sections are not changed. 

7 Conclusions 

1. Zee model reproduces rather naturally the pattern of neutrino masses and mixing 
which solves the atmospheric neutrino problem, supplies a desired HDM component 
in the Universe and gives the z/^ —T7 e oscillations in the range of sensitivity of existing 
experiments. 

2. The solar neutrino problem can be solved in extension of the model with an 
additional singlet fermion s, so that solar neutrinos undergo u e — > s conversion. The 
introduction of s does not destroy basic features of the Zee model. 

3. The data on oscillations of solar and atmospheric neutrinos as well as the 
cosmological mass scale fix all parameters of the Zee mass matrix. The scenario 
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implies in general inverse flavour hierarchy of the Zee boson couplings. There is a 
possibility of f efl ~ f^ T ^> f eT which may imply certain horizontal symmetry. 

4. The masses of the charged scalar bosons are of the order 100 - 500 GeV, and 
in certain cases at least one of the bosons can be as light as the lower kinematical 
bound. 

5. The scenario will be tested in forthcoming experiments: 

(i) the probability of — > u T oscillations is expected to be very small; discovery of 
these oscillations in CHORUS and NOMAD will rule out the scenario, (ii) The signal 
of v e — > v T oscillations may be in the region of sensitivity of these experiments, (hi) 
The confirmation of the LSND positive result will further testify for the suggested 
scenario, (iv) One may expect deviations from the SM predictions in ji — > v^eVf, 
decay, e~ — > e~ and z/ M e~ — > z/ M e~ scatterings, the violation of e — fi — r 
universality etc.. (v) The /i — > cy decay can be close to the present experimental 
upper bound, (vi) The life time of the neutrino radiative decay z/ 3 ( 2 ) — * ^17 is ex- 
pected to be 10 22 — 10 29 years. The decay of the relic neutrinos may have observable 
astrophysical consequences. 
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Figure Captions 



Fig.l The dependence of the oscillation probabilities P{y e — > v T ) at CHORUS 
and NOMAD (solid line), and Piy^ — > f e ) at LSND (dashed line) on sin6> for 
Am 2 2 = 10~ 2 eV 2 and m = 2.45 eV. 

Fig. 2 The ratios / MT / ' f eT (solid line) and fen/ fer (dashed line) as the functions 
of m for Piyp — > i/ e ) = 1.5 x 10~ 3 . 

Fig. 3 The lower bound for sin0 as the function of M\ for Piv^ — > z/ e ) = 
3 x lCT 3 (solid curve), 3 x 1CT 4 (long-dashed curve) and 3 x 1CT 5 (short-dashed curve). 
For other parameters we take (a) tan/5 = 59.3, m = 2.45 eV, (b) tan/3 = 2, 
m = 2.45 eV (c) tan/3 = 59.3, m = 10 eV. In all the cases M 2 = 300 GeV. 

Fig. 4 The dependence of the branching ratio of /i — > e^y on m for sin0 = 0.02 
(solid curve), 0.05 (dashed curve), and 0.1(dashed-dotted curve). The values of 
other parameters are fixed as tan (3 = 20, M x = 500 GeV, M 2 = 300 GeV and 
sin 2 29 = 2 x 10~ 3 . The short vertical lines indicate the lower bounds on m from 
the muon decay. The experimental upper bound on B(/i — > ej) is shown by the 
horizontal dashed line. 
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